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It has been argued that the quantum (conformal) trace anomaly could potentially provide us with 
a dynamical explanation of the cosmological constant problem. In this paper, however, we show by 
means of a semiclassical analysis that the trace anomaly does not affect the cosmological constant. 
We construct the effective action of the conformal anomaly for flat FLRW spacetimes consisting 
of local quadratic geometric curvature invariants. Counterterms are thus expected to influence the 
numerical value of the coefficients in the trace anomaly and we must therefore allow these parameters 
to vary. We calculate the evolution of the Hubble parameter in quasi de Sitter spacetime, where 
we restrict our Hubble parameter to vary slowly in time, and in FLRW spacetimes. We show 
dynamically that a Universe consisting of matter with a constant equation of state, a cosmological 
constant and the quantum trace anomaly evolves either to the classical de Sitter attractor or to 
a quantum trace anomaly driven one. When considering the trace anomaly truncated to quasi 
de Sitter spacetime, we find a region in parameter space where the quantum attractor destabilises. 
When considering the exact expression of the trace anomaly, a stability analysis shows that whenever 
the trace anomaly driven attractor is stable, the classical de Sitter attractor is unstable, and vice 
versa. Semiclassically, the trace anomaly does not affect the classical late time de Sitter attractor 
and hence it does not solve the cosmological constant problem. 

PACS numbers: 98.80.-k, 04.62. +v, 95.36,+x 



Recent observations have clearly indicated that the expansion of the Universe is accelerating. According to Einstein's 
general relativity, this can only be realised if the pressure of the dominant component of the current Universe is 
negative. These observations have triggered a renewed interest in the cosmological constant problem (for recent 
reviews, see e.g. [3, El El ) . What is usually referred to as the "old" cosmological constant problem can be phrased as 
follows: why is the measured (effective) cosmological constant extremely close to zero? 

One approach dealing with the cosmological constant problem is concerned with employing the effective field theory 
of gravity [3, . Lacking a full quantum theory of gravity, an effective field theory of gravity adopts the following point 
of view: in order to describe quantum phenomena at very large and cosmologically relevant distances, the precise 
physics at the shortest distance scales is irrelevant. In other words, the effective field theory of gravity is the low 
energy limit of quantum gravity. It combines classical general relativity with knowledge of quantum field theory in 
curved spacetimes @ . 

In order to describe these long distance effects accurately, one supplements the classical Einstein-Hilbert action with 
certain additional contributions. One of these additions is the trace anomaly or conformal anomaly which quantum 
field theories are known to exhibit @, 0, H, H, EH EH El, EH- If the classical action is invariant under conformal 
transformations of the metric, the resulting stress-energy tensor is traceless. As an explicit example, one can easily 
verify that the trace of a massless, conformally coupled scalar field vanishes. In quantum field theory the stress-tensor 
is promoted to an operator. A careful rcnormalisation procedure renders its expectation value (T 7 ^) finite. However, 
inevitably, the renormalisation procedure results in general in a non- vanishing trace of the renormaliscd stress-energy 
tensor. Classical conformal invariance cannot be preserved at the quantum level. Ever since its discovery, the trace 
anomaly has found many applications in various areas in physics (see e.g. [Hj]). 

An alternative approach to the backreaction problem of quantum fluctuations on the background spacetime deals 
with quantum fields whose spectrum is nearly flat [H, EH ES El, EH EH EH, El El, EI El|. Consequently, the 
spectrum in the infrared is not suppressed and is therefore expected to yield a strong backreaction on the background 
spacetime. Examples of such fields are the minimally coupled massless scalar and the graviton. 
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A. The Connection between the Cosmological Constant and the Trace Anomaly 

Some authors stated that the trace anomaly could have effects on dark energy and the cosmological constant 
problem (2|| [27[ , whereas it has been argued by other authors that the trace anomaly could potentially provide us 
with a dynamical explanation of the cosmological constant problem [28l, [29l , l30l . l3ll . [32 . l33j . Broadly speaking, the 
line of reasoning is as follows (for a more in-depth review, we refer to [33j|). The new, conformal degree of freedom is 
usually parametrised by: 

g liV {x) = e 2 °^g llv {x). (1) 

According to the authors of e.g. (33[, the trace anomaly cannot be generated from a local finite term in the action, 
but rather stems from a non-local effective action that generates the conformal anomaly by variation with respect 
to the metric [34] . It is this genuine non- locality of the trace anomaly, revealing a large distance effect of quantum 
physics, that is at the very foundation of its connection with the effective field theory of gravity. One then argues that 
the new conformal field should dynamically screen the cosmological constant, thus solving the cosmological constant 
problem. 

B. The Semiclassical Approach to the Cosmological Constant and the Trace Anomaly 

The proposal advocated in [33| is very interesting and should be investigated further. Before studying the effect of 
a new conformal degree of freedom (TT]) , we feel that firstly a proper complete analysis of the dynamics resulting from 
the effective action of the trace anomaly should be performed. This is what we pursue in this paper. 

According to the Cosmological Principle the Universe is homogeneous and isotropic on the largest and cosmologically 
relevant scales. The CMB measurements [HI constrain the inhomogeneities at order 10~ 4 ~ 1CU 5 . Moreover, the 
Universe appears to be spatially flat. Let us make the following observations. 

Firstly, the Cosmological Principle dictates the use of the conformally flat FLRW metric g^ u = a 2 (r))ri^. Hence, 
inhomogeneous fluctuations of the metric tensor and in particular of the conformal part of the metric tensor ([1]) are 
observed and expected to be small at the largest scales, comparable to the Hubble radius (also in the early Universe). 

Secondly, we are led to an essentially semiclassical analysis. The vacuum expectation value of the stress energy 
tensor resulting in the trace anomaly has been calculated semiclassically. In a semiclassical analysis quantum fluc- 
tuations backreact on the background spacetime. Phase transitions aside, quantum fluctuations naturally affect the 
homogeneous background homo gen ously. It is a well-known fact that quantum fluctuations can break certain sym- 
metries present in de Sitter [Til. Il6l ITtI Il8|. e.g. time translation invariance. However, we are not aware of quantum 
fluctuations breaking the homogeneity and isotropy of the background spacetime. 

Finally, if quantum fluctuations compensate for or screen the cosmological constant, we must have oc Q^u- Let 
us set: = 6{x)g liu . Stress energy conservation and metric compatibility immediately yield: V ^{x) = d fl 8(x) = 0. 
Hence we conclude that 9(x) must be a constant: 9{x) = 0q- Only homogeneous vacuum fluctuations can compensate 
the cosmological constant. Moreover, = B^g^ does not break any of the symmetries of a maximally symmetric 
spacetime 1 . Hence, this form cannot be used to study dynamical backreaction. 

The arguments above motivate a semiclassical approach to examining the connection between the cosmological 
constant and the trace anomaly. Note that we do not consider a new, conformal degree of freedom (fT]). Hence, we 
do certainly not exclude any possible effect this (inhomogeneous) conformal degree of freedom might have on the 
cosmological constant. However, it is plausible that in order to address the link between the cosmological constant 
and the trace anomaly, a semiclassical analysis suffices. 

C. The Modified Starobinsky Model 

Another application of the conformal anomaly can be found in what has become known as trace anomaly induced 
inflation: in the absence of a cosmological constant, the trace anomaly could provide us with an effective cosmological 
constant. Originally, Starobinsky [361 ] realised that quantum one loop contributions of massless fields can source a 
de Sitter stage. Subsequently, the theory of trace anomaly induced inflation received significant contributions from 
[37l . [H, HH El- If one includes a cosmological constant, the theory of anomaly induced inflation is plagued by 



1 Maximally symmetric spacetimes are de Sitter, anti de Sitter and Minkowski spacetime. 
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instabilities, which we will also come to address. The Modified Starobinsky Model as advocated by [HI, EH, IHj] , takes 
advantage of these instabilities to account for a graceful exit from inflation. It is argued that supersymmetry breaking 
changes the degrees of freedom such that it destabilises the quantum anomaly driven attractor and simultaneously 
stabilises the classical de Sitter attractor. 

Improving on e.g. (42l. we incorporate matter with a constant equation of state in the Einstein field equations. 
We argue that it is simply inconsistent not to include matter. Consider the following analogy: if we examine an empty 
Universe with a cosmological constant only, there is no dynamics and the (00) Einstein equation yields H 2 = A/3. 
Matter drives the dynamics and H 2 7^ A/3 can only be realised with pm 7^ 0. 

In the literature, if one solves the trace of the Einstein field equations in a Universe with a cosmological constant 
and trace anomaly, one solves, however, in reality for the dynamics in a Universe filled with radiation. For it is only 
radiation with equation of state w = 1/3 that does not contribute to the trace of the Einstein field equation T ra( j = 0. 
This point has not been included in other papers. We will consider matter with constant but otherwise arbitrary 
equation of state w > — 1, and not just (implicitly) radiation with w = 1/3. 

D. Outline 

In this paper we show that the cosmological constant problem cannot be solved by taking account of the trace 
anomaly alone. The outline of this paper is as follows. In section [TT] we recall the basics of the conformal anomaly 
and discuss how to study its effect on the evolution of the Universe by tracing the Einstein field equations. 

In section IIII1 we derive the conformal anomaly from an effective action in flat homogeneous FLRW spacetimes 
consisting of local quadratic geometric curvature invariants. Since one usually adds infinite counterterms to cancel the 
radiative one loop divergences, we do not see any reason why we should exclude adding a Gauss-Bonnet countertcrm 
to cancel the anomaly in flat FLRW spacetimes. Even though this term in the effective action is formally divergent, at 
the level of the equation of motion it yields a finite result. Hence, the coefficients multiplying the curvature invariants 
in the trace anomaly are not uniquely specified by the anomaly. The physical coefficient, i.e.: the parameter that can 
be measured, receives contributions both from the trace anomaly and from possible counterterms cancelling diver- 
gences from the underlying (and yet unknown) fundamental theory. This motivates varying the coupling parameters 
multiplying the curvature invariants in the anomaly. We can thus study all possible effects of the anomaly on the 
evolution of our Universe. 

In section HVl we study the evolution of a quasi de Sitter Universe in the presence of matter with constant equation 
of state, a cosmological constant and the trace anomaly. In quasi de Sitter spacetime we assume, loosely speaking, 
that the Hubble parameter is a slowly varying function of time. Effectively, we truncate the expression of the exact 
trace anomaly and discard higher order derivative contributions. 

In section [V] we generalise our analysis and study the evolution of an FLRW Universe again in the presence of 
matter with constant equation of state, a cosmological constant and the trace anomaly. We examine the exact trace 
anomaly and take all higher derivative contributions into account. As the dimensionality of the phase space increases, 
we must carefully perform a stability analysis of the late time asymptotes. 

II. TRACING THE EINSTEIN FIELD EQUATIONS AND THE TRACE ANOMALY 
A. The Conformal Anomaly in Four Dimensions in FLRW Spacetimes 

The trace anomaly or the conformal anomaly in four dimensions is in general curved spacetimes given by @, 0, [H| : 

T Q = (T^) = bF + b' (e - ^OR^j + b"UR, (2) 

where: 

E = *R fWK x*R,^ x = R^,xR^ X - 4iV^ + R 2 (3a) 
F = C^ KX C^ X = R^ vkX R^ kX - 2R ltv Rr + ii? 2 , (3b) 

where as usual R^ VK ,x is the Riemann curvature tensor, *R I1VK \ = s^apR" 13 kX /2 its dual, C^^a the Weyl tensor 
and R^y and R the Ricci tensor and scalar, respectively. Note that E is the Gauss-Bonnet invariant. The general 
expression for the trace anomaly can also contain additional contributions if the massless conformal field is coupled 
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to other long range gauge fields (see e.g. [||). Finally, the parameters b, b' and b" appearing in ([2]) are dimensionless 
quantities multiplied by h and are given by: 



6 
b' 



1 



{N S + 6N F + 12N V ) 
11 

360(4^) 2 {^ S + Y 



120(4tt) 2 
1 



i N s + —N F + 62N V 



(4a) 
(4b) 



where Ns, Np and Ny denote the number of fields of spin 0, 1/2 and 1 respectively (h — 1). It is important to note 
that b > whereas b' < in general. It turns out that the coefficient b" is regularisation dependent and is therefore 
not considered to be part of the true conformal anomaly. We take this into account and study the effect of b" on 
the stability of the solutions we are about to derive. For definiteness, we will assume that these parameters take 
their Standard Model values: Ns = 4, Np = 45/2 and Ny = 12. Note if we were to include right-handed neutrinos, 
Np = 24. One could also examine the numerical value of the coefficients (H} for the late time Universe. Today's 
massless particle is just the photon, hence Ny = 1, Ns = and Np = 0. 

Let us specialise to flat Friedmann-Lcmaitre-Robcrtson- Walker or FLRW spacctimes in which the metric is given 
by gap — diag (— 1, a 2 (t), a 2 (t), a 2 (i)) where a(t) is the scale factor of the Universe in cosmic time t. Recall that a 
conformal transformation leaves the Weyl tensor invariant. Hence, in FLRW spacetimcs F — 0. Given the FLRW 
metric one can easily verify that: 



R 2 



36 
12 
12 



H 2 + 4{HH 2 + H 4 ) 
H 2 +3(HH 2 + H 4 ) 
H 2 + 2{HH 2 + H 4 ) 



DR 



6 ii + 7HH + m 2 + Y1HH 2 



Hence, the exact expression for the trace anomaly in FLRW spacetimes in four dimensions reads 
T Q = 4b' |)/ + 7HH + AH 2 + 18HH 2 + 6iJ 4 | - 6b" j 



H + 7HH + 4H 2 + 12HH 2 } 



(5a) 
(5b) 
(5c) 
(5d) 

(0) 



To capture the leading order dynamics we work in quasi de Sitter spacctime and allow for a mildly time dependent 
Hubble parameter: 



H 
H 2 



constant <C 1 . 



(7) 



i.e.: we assume that e is both small and time independent. This would truncate the trace anomaly up to terms linear 
in H yielding: 



To = 246' 



'{ 



3HH 2 +iT 4 j - 72b" HH 2 



(8) 



We will examine both the exact form of the trace anomaly © and its truncated form ([H]). 

Truncating the expression for the trace anomaly is motivated by the following realisation. In general backgrounds 
we need a non-local effective action to generate the trace anomaly in the equation of motion. The non-locality at 
the level of the effective action corresponds to an expansion in derivatives at the level of the equation of motion. 
Generally, higher derivative contributions in an equation of motion have the tendency to destabilise a system unless 
the initial conditions arc highly fine-tuned. Formally, this is known as the theorem of Ostrogradsky and its relevance 
to Cosmology is outlined, for example, in [4oT ] . 

Note that when b" = 26'/3 the truncated version is exact. We discuss this further in section HVB[ Finally, note 
that although we have truncated equation ^ to obtain ([5]), equation ([5]) is still covariant. 



B. The Dynamics driven by the Trace Anomaly 

From the Einstein-Hilbert action: 

S = Seh + = 



16^G / d ^ x ^~ 9 ( R - 2A ) + / d 4 xV^gC M ■ 



(9) 
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where: 



Cm = -^d a <f>(x)d^(x)g^ - ^^(x) - V{<j>{x)) , (10) 



the Einstein field equations follow as usual as: 

Rp. v - ^Rg^u + ^g^ v = SirGT^ , (11) 

of which the trace can easily be verified to be: 

R - 4A = SttGT , (12) 

where T = T^. If one considers an empty Universe with a cosmological constant, the (00) Einstein field equation 
acts as a constraint equation for the Hubble parameter and one simply finds H 2 = A/3 as usual. However, for a 
non-empty Universe, the (00) Einstein field equation becomes a dynamical constraint. The Bianchi identity for the 
left-hand side of equation (fTTj) straightforwardly results in stress-energy conservation for the right-hand side: 

V%, = 0. (13) 

Because of stress-energy conservation, the (00) and the (ij) components of the Einstein field equations are not inde- 
pendent 2 . Therefore, any linear combination of the (00) and (ij) components of the Einstein field equations combined 
with stress-energy conservation suffice to describe the time evolution of the Hubble parameter. In particular, the 
trace equation (|12p and stress-energy conservation Q13p contain all relevant dynamics for H. 

Let us set <f)(x) = (j)ci{x) + f{ x ) f° r the quantum field in Sm and require that the classical field obeys the equation 
of motion. Note that the quantum perturbation ip(x) does not obey this equation of motion. We expand in terms of 
the quantum field and construct the effective action as usual: 

/( f & 1 f b 2 S 

XVexp i( / f M x) + - / ——-^——Lp(x)ip(y)+0((p 3 ) 
L \Jx 0(p c \{x) 2 J x y b0 c \{x)^ c \{y) 

= eaq>[iS M [&d+ir Q [<ki]] . (14) 

The first contribution to the effective action corresponds to the classical part of the action and Tq^i] is the con- 
tribution to the effective action taking account of the vacuum fluctuations. The stress-energy tensor now follows 
as: 

V = - 7=^ r ^ = - 7=f^ (5m[0ci] + rQ[0cl]) s T £ + T ^ ■ (15) 

Hence, there are both classical and quantum contributions to the full stress-energy tensor. Classically, from the 
equation of motion the scalar field obeys, we have: 

Vl£ = 0. (16) 

Hence, from (|13|) we derive: 

V% Q „ = . (17) 

Concluding, due to stress-energy conservation at the classical level and for the full stress-energy tensor, we have 
derived stress-energy conservation for the quantum contributions as well. 

Analogously to the classical stress-energy tensor, we can symbolically write: T' i u q = (— /Oq,Pq,Pq,Pq). Combining: 

T Q (t) = - PQ (t)+3p Q (t) (18a) 
pq(t) = -3H{p Q (t)+p Q (t)} , (18b) 



2 For example, stress-energy conservation combined with the (00) Einstein field equation straightforwardly yield the (ij) component of 
the Einstein field equations. 
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yields: 



dt 



[a 4 (t)p Q (t)]=-a 4 (t)H(t)T Q (t). 



We thus find (identical to |37j): 



PQ(t) = dra 4 (r)H(r)T Q (r) 



PQ(*) 



(T Q (t) + p Q (t)) 



(19) 

(20a) 
(20b) 



Although in general spacetimes it is not possible to perform this integral, in cosmologically relevant FLRW spacctimes 
we can [4y]. If we consider the exact form of the trace anomaly © in flat FLRW spacetimes, we can easily see that 
PGi(t) should be of the following form: 



p Q (t)=b' Cl HH + c 2 H 2 + c 3 HH 2 + c A H 4 



b" 



c 5 HH + c 6 H 2 + c 7 HH< 



(21) 



Upon inserting this ansatz into equation (|19[) and equating the contributions at each order, we immediately find: 



PQ 



2b' 



-2HH + H 2 - 6HH 2 - 3H 4 



3b" 



2HH - H 2 + 6HH 2 



(22) 



Because pq can be expressed in a local form, we would like to point out that this yields a local expression for the stress- 
energy tensor too. Although quantum fluctuations of the conformal field may still act non-locally on the background 
spacetime, at the classical level the trace anomaly affects the spacetime only locally. Again, when working in quasi 
de Sitter spacetime where e is both small and time independent, we can truncate this expression for the quantum 
density finding: 



p Q = -66' 2HH 2 + H 4 



18b" HH 2 . 



(23) 



Again, when b" = 2b' /3, this analysis becomes exact. 

In the next sections, we solve the Einstein field equations in FLRW spacetimes for a Universe in the presence of 
a) a non-zero cosmological constant, b) the trace anomaly as a contribution to the quantum stress-energy tensor and 
c) matter with constant equation of state pyi = wpM, where w > — 1. We thus consider non-tachyonic matter only, 
which docs not exclude [47], EH, where the effect of a finite period with w < — 1 is investigated. 

The trace anomaly enters the Einstein field equation naturally in the trace equation (|12p . We can thus study the 
effect of the trace anomaly on the evolution of a Universe with a cosmological constant and matter. Hence, the trace 
equation of the Einstein field equations is given by: 



R — 4A = —8irG {Tq + Tm} 



(24) 



where 



T M = ~Pm + 3p M = Pm (3w - 1) . (25) 

Now, we can employ the (00) Einstein field equation from (jlip to express pm in terms of pq yielding: 

9(1 + u)H 2 (t) + 6H(t) - 3(1 + w)A = -8ttG [Tq + (1 - 3lo) Pq ] . (26) 

The above equation governs the dynamics for the Hubble parameter that we will solve in various interesting cases. For 
the anomalous trace we can either take the exact expression ([5]) containing higher derivatives of the Hubble parameter 
or its truncated version ([8]). Likewise, for the quantum density we can either insert the full expression (|22[) or the 
truncated one (|23p . The reason for truncating the expression for the quantum trace and density as outlined above is 
the realisation that higher derivative contributions in an equation of motion generally have the tendency to destabilise 
a system. 

In the literature (see e.g. [HI), one only has considered an empty Universe with a cosmological term and trace 
anomaly. In reality, one solved for a radiation dominated Universe, for only radiation (lo = 1/3) does not contribute 
to the trace of the Einstein field equations. We incorporate matter with constant equation of state parameter ui. 
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Independently on whether one truncates the expressions for the anomalous trace or quantum density, one can solve 
for the asymptotes yielding the late time behaviour. Setting all time derivatives of the Hubble parameter equal to 
zero, one can easily solve for the two late time constants: 

\ H ° ) = 32^ ' (27) 

Here, Hq turns out to be the classical de Sitter attractor, whereas Hq is a new, quantum anomaly driven attractor. 
This result for the late time constants is identical to the case where matter is absent [42j |. We can write the above 
expression in a somewhat more convenient form by defining the dimensionless parameter A: 

A = — , (28) 

that sets the scale for the cosmological constant A. In the current epoch, A is extremely small, which allows us to 
expand (|2"7|) finding 3 : 

#o C = y|[l-87r6'A] (29a) 



< = Vl6^-I- (2%) 

In the absence of a cosmological constant, the trace anomaly can thus provide us with an inflationary scenario which 
has already been appreciated by [H, H3, OH, H3| ■ Finally, note these asymptotes are independent of b" . 



III. THE EFFECTIVE ACTION GENERATING THE TRACE ANOMALY 

The authors of [33| are correct in saying that the conformal anomaly cannot be generated by a finite effective action 
built out of local quadratic geometric curvature invariants only. We show that the trace anomaly can be generated 
by an infinite effective action in flat FLRW spacetimes, that consists only of local quadratic geometric curvature 
invariants. Although infinite at the level of the effective action, we generate a finite on shell contribution. We first 
write the trace anomaly in conformal time dt = a(rj)dr] such that the full conformal anomaly ([6]) reads: 

,,(a"{a'\ 2 fa'V\ ,,,,fa"" a'" a' fa"\ 2 a" { a' \ 2 \ , x 

r « = 245 U) - b) J - 6 v - 26 ^ [ur ~ - 3 U) + 6 ^ UJ J • (30) 

Here, dashes denote conformal time derivatives. In general, the trace of a stress-energy tensor can be written as: 

r = — H-^"-^- = — L_ 6 A Of) 

Here, V is the (spatial) volume. The correct effective action T an that generates the trace anomaly in spatially flat 
FLRW spacetimes is given by: 

r an = J d D x^—g {/3 D E - 12 (b" - 2b' /i) R 2 } , (32) 



where 



/3d = 5 '0D^4) (33) 



3 Note the nomenclature in the literature is somewhat misleading. Rather than calling H29al the classical de Sitter attractor, it would 
be more natural to denote it with the quantum corrected classical attractor. Hence, the quantum attractor I l29bl l should preferably be 
denoted by anomaly driven attractor or Planck scale attractor. We will nevertheless adopt the nomenclature existing in the literature. 
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Here D is the dimension of the spacetimc and @q is an undetermined finite (and physically irrelevant) constant. Only 
after variation we let D — > 4. The numerical factors in are chosen in accordance with the trace of the Einstein 
field equations, i.e.: the trace of the variation: 

„u J_ S = g^ [S EH + r an ] = , (34) 

indeed yields: 

R - 4A = -BttGT . (35) 
The effective action (|32|) merits some clarifying remarks. Firstly, the variation of R 2 straightforwardly yields the OR 



term in the conformal anomaly. This term is not unique for one could have equally well taken R^ U R^ U into account 4 . 
Secondly, note that the local effective action depends solely on the scale factor. If one were to rewrite this effective 
action covariantly in terms of the full metric, the effective action would become non-local (8l. [3 3 . also see 0, [5(3] ■ It is 
well-known that the effective action of the conformal factor can be written in a local form [46|, [51(. However, we have 
been able to rewrite this expression in terms of the Gauss-Bonnet invariant multiplying an infinite constant. Finally, 
note that the coefficient (|33[) diverges in four dimensions. However, the Gauss-Bonnet invariant in D dimensions 
reads: 

which can be determined from the D dimensional generalisations of equation Note: 

D{D-A)E. (37) 



Clearly, only in four dimensions the Gauss-Bonnet term corresponds to a total derivative that can be neglected. 
Hence, we choose the divergent coefficient /3d such that the factor (D — 4) cancels yielding a non-vanishing and finite 
contribution at the level of the equation of motion when D — > 4. The contribution from /3q is identically zero when 
we let D — > 4. Although divergent in the action, we have shown that at the level of the equation of motion the term 
involving (3d yields a finite contribution. All physical measurements arc performed on shell, hence we cannot exclude 
a counterterm of the form (3dE in the effective action. 

This procedure can be debated. In the literature one usually prefers a finite effective action. We abandon this 
assumption and require a finite equation of motion only. We cannot think of a physical measurement that distinguishes 
the two approaches and therefore we argue that one should consider and examine all possible effects the countertcrms 
might have on the trace anomaly. 

In homogeneous cosmology, R 2 and E are the only local quadratic geometric curvature invariants 5 . Divergences 
up to one loop in perturbative quantum gravity can be cancelled only by counterterms of the form R 2 and E. The 
physical coefficients multiplying R 2 and E in the rcnormaliscd one loop effective action receive contributions which 
one can write as: 

aphys-R 2 = [a anom {m) + a ct ] R 2 (38a) 
P P hy S E = [(3 &nom (m)+f3 ct ]E, (38b) 

where a anom (m) is the mass dependent finite contribution from the trace anomaly, where m = {rrii} denotes the 
mass of the particles i. It is in principle uniquely defined by the requirement that a an0 m — » as m — * oo. Similarly, 
/?anom(w) is the (infinite) anomalous contribution determined by the requirement that /3 a nom — ► as m — > oo. These 
requirements in principle fix /3 an om(?Ti) and a a nom(™) uniquely. While there is agreement in the literature on the 
value of a a nom(™) (also used in this paper), disagreement exists on the value of f3 anom (m). Therefore we leave it 
unspecified. The contributions a c t and /3 c t correspond to the parts of the counterterms that remain when eventual 
one loop divergences are cancelled 6 . It is however only the sum of these terms, yielding a p h ys and /3 p hys, that is 



4 We can easily split in an B? contribution and a Gauss-Bonnet term, that yields a vanishing surface term in four dimensions. 

5 Because F = 0, we can express R I1 , VK \R^ VK ' > ' in terms of R llv R> iv and R 2 . Hence, R? and E are linearly independent. Note furthermore 
that if one investigates inhomogeneities in the Universe, one measures (statistical) correlation functions. They are translation invariant 
as a consequence of the symmetry of the vacuum state and therefore respect the symmetry of the background spacetime. 

6 For a calculation involving anomaly calculations around Minkowski spacetime, we refer to [52I, |53||. For the calculation of the o anom ('7x) 
function for a scalar field in de Sitter spacetime, we refer to [f|. 
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physical, i.e. measurable [2J, |25|. Since these parameters have not been measured, we cannot simply assume that the 
coefficients in the trace anomaly are just given by ([!]). We should allow these coefficients to vary in order to examine 
the full effect of the trace anomaly on the evolution of our FLRW Universe. This is precisely what we pursue in the 
following sections. 



IV. THE TRACE ANOMALY IN QUASI DE SITTER SPACETIME 

In this section we work in quasi de Sitter spacetime, where we treat e as a small and time independent constant 
which allows us to neglect higher order derivative contributions. The reason for discarding higher order derivative 
contributions is that they tend to destabilise a dynamical system, formally known as the theorem of Ostrogradsky 
0. 

We have to distinguish two cases separately. In the spirit of [33j , the numerical value of the parameter b" occurring 
in the trace anomaly is not fixed because it is a rcgularisation scheme dependent parameter. Generally, however, we 
cannot exclude the presence of this term and we therefore allow it to take different values. First, we allow for an 
unrestricted value of b" and secondly we set b" = 26'/3. This case is particularly interesting as this value of b" sets 
the total coefficient multiplying the OR contribution in the trace anomaly to zero. 



A. Case I: unrestricted value of b' 



We thus insert the truncated expression for the trace anomaly ([5]) and the quantum density (|23p into the Einstein 
field equation (|2l))) . This yields: 



9(1 + uj)H 2 (t) + 6H(t) - 3(1 + uj)A = -8nG {126'(5 + 3w) - 546"(1 + u)} HH 2 + 186'(1 + uj)H 1 
This differential equation can be solved exactly. Separation of variables yields: 
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where a is conveniently defined as: 



a = 8irG {46'(5 + 3w) - 186"(1 + u)} 



(39) 
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One can perform the integral in terms of logarithms where one has to take the signs of the occurring parameters 
carefully into account. The integral above gives: 
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(42) 
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Note the asymptotes of this analytic solution coincide with the asymptotes obtained earlier in (f2"9"]) as expected. 

In figure [l] we numerically calculate the dynamics of the Hubble parameter for various initial conditions. The two 
asymptotes divide this graph into three distinct regions that are not connected for finite time evolution. The region 
bounded by the two asymptotes contains initial conditions for H(t) such that H{t) grows for late times towards Hq 
and initial conditions such that H(t) asymptotes to the de Sitter attractor Hq. In figures [5] and [3J we examine 
whether our approximation that e is both small and a constant is valid for calculating the dynamics. In figure [2] one 
can clearly see that e 1 for late times. If e = 0, we should have e/(He) <C 1, an assumption that is violated as 
depicted in figure[3J However, already for a classical cosmological constant dominated Universe with matter, a similar 
violation occurs. 

Furthermore, note the existence of a branching point, an initial condition for H (t) such that for H(0) > Hbp the 
Hubble parameter asymptotes to the quantum attractor Hq and for H (0) < Hbp the Hubble parameter decreases to 
the classical de Sitter attractor Hq. When rewriting equation of motion (|3"9"|) in terms of e and noting that e should 
diverge exactly at the branching point, one finds: 
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Figure 1: Dynamics of the Hubble parameter in quasi de 
Sitter spacetime in the presence of a non-zero cosmolog- 
ical constant, the trace anomaly, and matter (u) = 0). 
Depending on the initial conditions, the Hubble parame- 
ter evolves to either the classical de Sitter or the quan- 
tum anomaly driven attractor. We have used X — 1/50, 
b' = -0.015 (Standard Model value) and b" = 0. 
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Figure 2: Validity of the assumption e <C 1. For the 
various initial conditions in figure QJ one can clearly see 
that this approximation is well justified. 
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Figure 3: Validity of the assumption e = 0. For the 
various initial conditions in figure QJ we have calculated 
e/[He). Clearly, this condition is violated at all times. 



Figure 4: Instability of the quantum anomaly driven at- 
tractor in quasi de Sitter spacetimes. We have used 
X = 1/50, w = 1/3, b" = 7fe'/6 < 56'/6. 



Wc probe the dependence on scale by changing the numerical value attached to A in equation (|28|) . If we decrease 
A, then also A decreases which results in a smaller Hq. Also, it turns out that both of the asymptotes are already 
reached much faster. This improves the validity of the assumption e <C 1, whereas the assumption e/(He) <§C 1 is still 
seriously violated at all times. 

Let us study the analytical solution more closely. In particular, it is interesting to derive the high and low 
energy limits of this solution. Of course, one naively expects in the high or low energy limit to flow towards the 
quantum or classical attractor, respectively. However, the analysis turns out to be somewhat more subtle. We will 
show that under a certain condition, the quantum anomaly driven attractor becomes unstable. Although solution 
(|42|) looks complicated at first glance, it simplifies when defining: 
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Wc define the initial conditions at t' as: 



(-1 
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H(f) 
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H{t>) 



(44d) 
(44e) 



Note that fi > 1, gencrically. With these definitions, equation 



n(t-t') = -a 

In the high energy limit, we set: 

S(t) = 

such that S(t) <C 1. Equation (j4"5)) thus modifies to: 
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(47) 



We can expand the second logarithm making use of A <C 1. The leading order contribution (in 8{t)) is given by the 
denominator in the logarithm, because this term diverges as 6(t) approaches zero. We can thus exponentiate the 
equation and solve for the Hubble parameter: 



\H>(t)-H^\ = 2iJ^exp 



where the time shift Ai> is given by: 



H c 

rtAt> =B[c 2 - 2-^ 



(t-t' + At>) 



(48) 



(49) 



As ft > 0, this solution converges whenever A < 0. This provides a stability condition on b" in terms of b' and A. 
The quantum anomaly driven attractor is stable, whenever the following inequality is satisfied: 



,.. 2,,4 + 3w + 167rA6'(5 + 3w) 

b > -b 

9 (1 + w)(1 + 16ttA6') 



(50) 



In figure 01 we have numerically calculated the evolution of the Hubble parameter in a radiation dominated Universe 
when this inequality is not satisfied. We used b" = 7b' /6 < 56'/6. For initial conditions above the attractor, the Hubble 
parameter increases to even higher energies, whereas for initial conditions below the quantum attractor, the Hubble 
parameter evolves towards the classical attractor. Hence even in quasi de Sitter spacctimes, physically questionable 
solutions occur. However, note that when b" = 26'/ '3, the specific case under consideration in subsection IIVB) the 
above inequality is satisfied. 

When the above inequality is satisfied, the Hubble parameter in the high energy limit decays exponentially towards 
the quantum anomaly driven attractor, where some "frequency dependence" through Q/A and a time shift At > can 
be recognised. The time shift can without observational consequences be absorbed in the initial time t' . 

The low energy limit reveals less surprising behaviour. Here, we set: 



S(t) = 



and S(t) <C 1. Again we use A 1 in order to capture the leading order dynamics. This yields: 
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Figure 5: Phase space flow in quasi de Sitter spacetime. 
The phase space consists of two lines. The flow, indi- 
cated by the arrows, is towards the classical or quantum 
attractor represented by two dots. In this regime, both at- 
tractors are stable. The vertical dashed line indicates the 
branching point and the dashed-dotted line the classical 
evolution. We have used u = 0, A = 1/50, b' — —0.015 
(Standard Model value) and b" — 0. 



Figure 6: Phase space flow in quasi de Sitter spacetime 
for b" — 2b'/3 such that the OR does not contribute in 
the trace anomaly. Qualitatively, the dynamics does not 
change when compared to figure [5[ Again, both attractors 
are stable. Apart from b" , the value of the parameters are 
identical to figure [5l 



Because B > 0, this solution converges. In the low energy limit the time shift At< is slightly different as compared 
to (HJ: 

flAt < =Bc 2 -A \ci - . (53) 

Finally, we examine the phase space flow in quasi de Sitter spacetime. In figure [5J we show a parametric plot of 
H versus e. The phase space basically consists of two lines. The phase space flow is towards either the classical or 
the quantum anomaly driven attractor as indicated by the arrows. Note that we have chosen both attractors to be 
stable. Furthermore, we also include the branching point (|13"|) and the classical evolution, that is, the evolution of a 
Universe with b" = b' = 0. As expected, the flow is towards the classical attractor in this case. Although the analysis 
performed above is for generic values of 6", we set it to zero in figure [3] and b' takes its Standard Model value. 



B. Case II: b" = 2b'/3 

As indicated earlier, we must consider the case when b" = 2b' /Z separately because in this particular case the total 
coefficient in front of the OR contribution to the trace anomaly vanishes. All higher derivative contributions precisely 
cancel and also the H 2 contribution happens to cancel, such that we find ourselves immediately situated in quasi de 
Sitter spacetime. Albeit a simple case, we do take the full trace anomaly into account. 

The analytic solution obtained in (|42[) still applies and moreover, it becomes exact. The branching point is still 
given by equation for which we just have to insert b" = 2b' /Z. Clearly, in figure [S] one can see that qualitatively 
the dynamics has not changed compared to figure [5j The branching point has shifted somewhat to the right, and the 
way in which the Hubble parameter approaches its two late time asymptotes differs. However, the important features 
of figure [5j i.e.: two stable attractors, the occurrence of a branching point and the shape and dimension of the phase 
space, do not change. 



V. THE TRACE ANOMALY IN FLRW SPACETIMES 



We turn our attention to solving the full trace equation (f2"o]l . where we truncate the expression neither for the 
anomalous trace © nor for the quantum density (|22[) . Obviously, we cannot solve this equation analytically, for it 
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Figure 7: Dynamics of the Hubble parameter taking the 
full trace anomaly into account. We took b =0 such that 
b" — 2b' / 3 > yielding an unstable classical attractor. We 
have used ui = 0, A = 1/50, b' — —0.015 (Standard Model 
value). 




Figure 8: Parametric phase space plot for figure [7] in 
which the classical attractor is unstable. We have indi- 
cated the quantum anomaly driven attractor as a small 
black sphere. 
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Figure 9: Dynamics of the Hubble parameter taking the 
full trace anomaly into account. We took b" = b' such that 
b" — 2b'/3 < yielding a stable classical attractor. We 
have used w = 0, A = 1/50, b' = -0.015 (Standard Model 
value). Clearly, the classical attractor is under-damped, 
resulting in various oscillations around Hq . 
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Figure 10: Parametric phase space plot for figure^ in 
which the classical attractor is stable. Because of under- 
damping, one spirals towards the classical attractor. For 
clarity, we have only included the flow for two initial con- 
ditions H(0) = andH(0) = 6y/W/A. The twister 
like structure is clearly visible. Qualitatively, the phase 
space flow resulting from other initial conditions is iden- 
tical. 



contains all higher order derivative contributions, which forces us to rely on numerical methods. 

Firstly, we note that the asymptotes (|29p do not change by including the higher derivative contributions. Keeping 
Ostrogradsky's theorem in mind, we expect to incur all kinds of issues related to the stability of our system and 
asymptotes in particular. It is therefore essential to perform a stability analysis for small perturbations SH(t) around 
both of the asymptotes. We insert: 

H(t) = H^ A + SH(t) , (54) 

in equation (|26|) . where Hq' A can either denote the classical or the quantum attractor. For the small perturbations 
around these attractors we make the ansatz SH(t) = c cxp[££] . Wc linearise the trace equation finding the characteristic 
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equation from which we determine the eigenvalues £ of our system: 

2 (a» - 3i//2) £ 3 + 6 (// - 3i//2) (2 + w)if C ' A £ 2 + ({6^(5 + 3lu) - 27v(l + uj)} (# aA ) 2 - 3^ £ 

+9 (V (^o C ' A ) 2 " lj (1 + ^)#o C ' A = , (55) 

where // = — 87rG6' and Z/> = SirGb". Remarkably, the solutions of this third order equation are simple: 

= -3iJ c - A (l+w) (56a) 

orrC,A 

e (2) = -— | 1" \/A (56b) 

= -£i|_-VA, (56c) 

where: 

H^ A Y (IO^ + qA . (57) 



4(2/i-3i/) 



Clearly, when Re£W < for i = 1,2,3, the corresponding attractor is stable. Since we only consider non-tachyonic 
matter, eigenvalue (|56ap is negative. However, a finite period in which w < — 1 is not excluded (see e.g.: [47l |48|V If 
one were to consider other equations of state than the simple linear one pm = wPm, this statement might no longer 
hold [54l | . Hence only for £( 2 ) when A > 0, we could encounter a potential instability. Surprisingly, the stability 
analysis does not depend on the equation of state lu because lu enters only through equation (|56a[) . The condition for 
instability thus reads: 

A>pf^f. (58) 

We can rewrite this equation to find: 

4 + 8ttG (i? C ' A ) 2 (106' + 96") ^ 72ttG (6" - 26'/3) (i? C ' A ) 2 • ( 59 ) 

In the expression above, we should read the inequality > or < whenever b" — 26'/3 > or b" — 26'/3 < 0, respectively. 
We can now insert either the classical or quantum asymptotes previously derived in equation (|29p and verify which of 
the two above inequalities is satisfied. Upon inserting the expression for the classical attractor, equation (|59[) yields: 

1 + 32ttA6' ^ , (60) 

where A = GA/3 <C 1 as before. Only the first inequality > will be satisfied. Hence we conclude that the classical 
attractor is unstable if b" — 2b'/ 3 > 0. The converse will be true if b" — 2b'/3 < 0. Likewise, equation (|59[) for the 
quantum attractor after some algebra reads: 



32ttA , 

1-— tffco, 



(61) 



Concluding, when using A< 1, we unambiguously find: 



tp ,// r.7/ /r. n ,1 f Classical attractor unstable . 
11 b — 2b 3 > 0, then < „ , , , (62a) 

' I Quantum attractor stable v ' 

„ ,// /n „ f Classical attractor stable s 

If b" - 2673 < 0, then { „ , , , , , , . (62b) 

' I Quantum attractor unstable v ' 

Let us first of all recall that it is precisely the combination b" — 2b'/ 3 that multiplies the OR contribution in the 
trace anomaly. This calculation thus proves the statements about stability made in e.g. [42j using the Routh-Hurwitz 
method. Our proof is more general because we include a constant but otherwise arbitrary equation of state parameter 
lu > — 1. Moreover, while the Routh-Hurwitz method can only guarantee stability of a solution (when certain 
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determinants are all strictly positive) , it does not tell anything about instability [42|, [55[ . Furthermore appreciate 
that the singular point in this analysis, b" — 2b'/ 3 = 0, or equivalently 2/J./3 — v = 0, immediately directs us to the 
quasi de Sitter spacetime analysis performed in section HVB( where all higher derivative contributions precisely cancel, 
rendering both attractors stable. 

Let us compare figures [7] and [9] In the former figure, we used b" = such that b" — 2b' / 3 > 0, yielding an unstable 
classical attractor. However, if H(0) < Hq the quantum anomaly driven asymptote is not an attractor and the Hubble 
parameter runs away to negative infinity. In the latter figure, we set b" = b' such that b" — 2b' / 3 < which gives us 
a stable classical attractor. Likewise, for initial conditions H (0) > Hq the de Sitter solution is not an attractor and 
the Hubble parameter rapidly blows up to positive infinity. 

In figure [HI we can observe another interesting phenomenon. In this case, the classical attractor is under-damped, 
resulting in decaying oscillations around the de Sitter attractor. In figure [7] these oscillations are not always present. 
The eigenvalues (|56p develop an imaginary contribution resulting in oscillatory behaviour whenever: 



We thus find: 



A < . (63) 



8ttG (#o C ' A ) 2 (10b' + 96") ^ . (64) 



The inequality < or > holds whenever b" — 2b' / 3 > or b" — 2b'/ 3 < applies, respectively. Again, we verify which 
of the two inequalities is actually satisfied. To study oscillations around a stable classical attractor, we should take 
the > inequality (there are no oscillations around an unstable attractor). We thus find: 

1 + 2ttA(106' + 96") > 0. (65) 

Clearly, this inequality is always satisfied because A <C 1. We thus conclude that whenever the classical de Sitter 
attractor is stable, oscillations occur. Furthermore, we can insert the quantum anomaly driven attractor in equation 
(|64[) . Now, we should use the < inequality in order to study oscillatory behaviour around the (stable) quantum 
attractor. This yields: 

-l-~-87rA(106' + 9b") <0. (66) 
Oscillatory behaviour around the quantum attractor thus occurs when: 

In figure [TTJ we depicted the parameter space (the b' versus b" plane) resulting in oscillatory behaviour around the 
quantum anomaly driven attractor. First of all, we should have b" — 2b' / 3 > yielding a stable quantum attractor. 
Of course, one cannot have oscillations around an unstable attractor. Secondly, note that b' < because of equation 
()4b[) . Finally, when the newly derived inequality (|67|) is satisfied, oscillations occur. These considerations divide the 
phase space into three regions: a region where oscillatory behaviour occurs, another region which results in critically 
or over-damped behaviour and a part of phase space that is forbidden, as shown in figure [TT] 

Let us now return to discussing the classical attractor that shows its oscillatory behaviour manifestly. We calculate 
the frequency of oscillations around this attractor. Taking the square root of (|57[) and extracting an i we find the 
(quantum corrected) frequency for oscillations around the classical attractor: 




27rA(10b' + 9b") 
BkG(W/3-V>) ■ (68) 

Note that this frequency is independent on the equation of state parameter to. The frequency of oscillations around 
the quantum attractor can be found analogously. 

Let us analyse the phase space in the case of a stable quantum and classical attractor subsequently. In figure 
we visualise the phase space flow for the former case parametrically in the H(t), e(t) and e(t) directions 7 . The small 



7 Note we are not able to include the fourth dimension of the phase space, e(t). However, also e(t) rapidly approaches zero as time elapses. 
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Figure 11: Parameter space for oscillations around the quan- 
tum anomaly driven attractor. The gray region (light shaded) is 
excluded either because b' > or because the quantum attractor 
has become unstable. The red region (dark shaded) in parameter 
space shows oscillatory behaviour, whereas the white region is 
either critically damped or over-damped. We used X = 1/50. 



black sphere denotes the quantum anomaly driven attractor. For the latter case, we include in figure [TU] the phase 
space flow for just two initial conditions for H(t) for clarity. The under-damped oscillatory behaviour results in the 
twister like structure visible in the H(t), e(t) and directions. 

Finally, we would like to point out that including all higher derivative contributions, which thus corresponds to 
solving the trace equation exactly, modifies the dynamics of the Hubble parameter significantly. Attractors that 
were stable in the absence of higher derivatives under certain conditions destabilise. The reason for this, clearly, is 
attributable to the presence of the OR term in the trace anomaly, generating these higher derivative contributions. 
We do not know whether or not incorporating the higher derivatives is a sensible thing to do. Usually higher order 
derivatives tend to destabilise a system signifying that some particular solutions are not physical. Therefore, in the 
spirit of Ostrogradsky's theorem, one can question whether the analysis where higher derivative contributions are 
discarded is correct, or the analysis taking the full trace anomaly into account. 



The trace of the Einstein field equations in cosmologically relevant spacetimes together with stress-energy conserva- 
tion completely captures the dynamics of the Hubble parameter. We have derived the trace anomaly from an effective 
action in spatially flat FLRW spacetimes. It consists of the local quadratic geometric curvature invariants R 2 and 
the Gauss-Bonnet term E. Because of counterterms that are supposed to cancel divergences of the as yet unknown 
underlying fundamental theory, we expect the coefficients in the trace anomaly to change. The physical value of 
each of these coefficients receives contributions both from the anomalous trace and from these counterterms. Because 
we do not know the physical value these parameters will take, we must allow them to vary in order to examine all 
possibilities. 

We have studied the dynamics of the Hubble parameter both in quasi de Sitter and in FLRW spacetimes including 
matter, a cosmological term and the trace anomaly. In quasi de Sitter spacetime, where we restrict the Hubble 
parameter to vary slowly in time, we find that for various initial conditions H(t) asymptotes either to the classical de 
Sitter attractor, or to a quantum anomaly driven attractor. We find a region in parameter space where the quantum 
attractor destabilises. Otherwise, both attractors are stable. 

In FLRW spacetimes we include all higher derivative contributions in the trace anomaly. We perform a stability 
analysis for small perturbations around the two asymptotes. For b" — 26'/3 > 0, the quantum attractor is stable and 
the classical de Sitter attractor is unstable. On the contrary, for b" — 2b' /i < 0, the quantum attractor is unstable 
and the de Sitter attractor becomes stable. The singular point in this analysis, b" — 26'/3 = 0, immediately directs 
us to quasi de Sitter spacetime in which the dynamics is much simpler. In this case, both attractors are stable. The 
classical de Sitter attractor always shows under-damped oscillatory behaviour and we calculate the frequency of these 
oscillations. We analyse the phase space of the quantum attractor and conclude there is some region in parameter 
space for which oscillations occur. 



VI. CONCLUSION 
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There is no dynamical effect that influences the effective value of the cosmological constant, i.e.: the classical de 
Sitter attractor. Based on our semiclassical analysis we thus conclude that the trace anomaly does not solve the 
cosmological constant problem. 

We have studied both the truncated and the exact expression of the trace anomaly in flat FLRW spacetimes. We 
do not know which of the two approaches is correct. Keeping Ostrogradsky's theorem in mind, higher derivative 
contributions usually have the tendency to destabilise a dynamical system. Discarding these higher derivatives and 
studying the trace anomaly in quasi de Sitter spacctime would thus seem plausible. 

Finally, one could wonder whether the quantum anomaly driven attractor is physical. The quantum attractor is 
of the order of the Planck mass M p \, so only when matter in the early Universe is sufficiently dense, H ~ 0(M p i). 
We then expect to evolve towards the quantum attractor. However, at these early times we also expect perturbative 
general relativity to break down. Hence, this attractor might even not be there or it may be seriously affected by 
quantum fluctuations. Quantum fluctuations present at that epoch might even induce tunnelling towards the regime 
where H (t) asymptotes to the classical attractor. 
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